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In density functional theory (DFT), the exchange-correlation functional can be exactly expressed 
by the adiabatic connection integral. It has been noticed that as A — > oo, the A~ x term in the 
expansion of W(\) vanishes. We provide a simple but rigorous derivation to this exact condition in 
this work. We propose a simple parametric form for the integrand, satisfying this condition, and 
show that it is highly accurate for weakly-correlated two-electron systems. 



In density functional theory (DFT) [I], the exchange- 
correlation functional E xc [n] is exactly expressed by the 
adiabatic connection [2] [3] formula: 



E xc [n] = I dXW[n](X), 
o 



(1) 



where A is a coupling constant that connects the Kohn- 
Sham system (A = 0) to the true system (A = 1), while 
keeping the density n(r) fixed. The integrand, W(X), 
contains only potential contributions to E xc . The shape 
of Ty(A) has been much studied in DFT [4]. For example, 
the success of hybrid functionals that mix some fraction 
of exact exchange with a generalized gradient approxi- 
mation (GGA) can be understood this way [5]. There is 
ongoing research to use the low density (A — > oo) limit as 
information in construction of accurate models of VF(A) 
[3 [TJ H] • Recently, the adiabatic connection formula has 
been used directly in functional construction [9]. 

The expansion of W(A) in the high-density (weak cou- 
pling) limit for finite systems is known to be [7J: 



W(X) = W + WqX 



■ as A — > 0, 



(2) 



where Wg = 2E^ L2 , with E^ L2 the second-order coef- 
ficient in Gorling-Levy perturbation theory [6] I10[ lllj . 
The expansion in the low-density (strongly correlated) 
limit is believed to be [7J [T5] : 

W{X)^W 00 + W' 00 X- 1/2 + A^oo, (3) 

where is defined as the coefficient of A -1 / 2 in the 
expansion above, and can be calculated from the 
strictly correlated electron (SCE) limit [13] ■ In addition 
to these expansions, by definition the exact VF[n](A) is 
known to satisfy the following scaling property [7 : 



W[n](X) = XWi[n 1/x ], 



(4) 



where ni/\(r) is the scaled density, defined by n 7 (r) = 
7 3 n(7r), < 7 < 00. In the equations above, one 
can show that Wq = E x , the exchange energy, and that 
Woo is finite [6]. The dependence on A -1 / 2 in the low- 
density limit is because correlation dominates here, and 
the Thomas-Fermi screening length is proportional to 



A 



-1/2 



In practical DFT calculations, W{X) must be approx- 
imated. However, any approximate VF(A) should satisfy 



several exact conditions, such as Eqs. ^ and Q. 
In the erratum to Ref. [7] , Seidl et al. concluded that for 
the ISI model (see below) , the spurious A 2 In A term in 
E c [n\] is due to the A -1 term in the expansion of VF(A) as 
A — > 00 [Eq. ([3j]. In a recent work [T2], this was proved 
rigorously, but only by calculating zero-point oscillations 
about the strictly-correlated limit. In this paper, we pro- 
vide a simple derivation and how this exact constraint 
affects approximate functionals. Throughout this paper, 
we use atomic units (e 2 = H = /x = 1) everywhere, i.e. all 
energies are in Hartrees and all distances in Bohr radii. 

Any A-dependence can always be expressed in terms of 
density scaling. Using the fundamental relation of Levy- 
Perdew equation |14j . one finds: 



iy W (A)^*M-7 2 ^(^ 



(5) 



and it is generally believed for nondegenerate Kohn-Sham 
systems |15] that E c [n 7 ] has the following expansion in 
the low density limit (7 — > 0) : 



E c [n 7 ] = 7 (SoM + 1 1/2 BM + 1 B 2 [n] 



(G) 



where the Bk[n]'s (k = 0,1,2 ■■■) are scale-invariant 
functionals. Substituting into Eq. ([5]), we find the ex- 
pansion of W(X) for large A: 

W(X) = E x [n} + B [n} + ^X- 1 / 2 B 1 [n}-^X-^ 2 B 3 [n} + - ■ ■ 

i.e. the A 1 term is missing, and T'F(A) is independent of 
B 2 [n]. 

Now we survey approximations to W(X) and see 
whether they have the correct low-density expansion [Eq. 
(firl]. There are several kinds of approximations, the most 
famous being the ISI (interaction-strength interpolation) 
model by Seidl et al [5J : 



VF 1M [n](AH W QO [n] + 



X[n] 



yJl + Y[n]X + Z[n] 



(8) 



where X = xy 2 /z 2 , Y = xX/z 2 , Z = X/z - 1, with 
x = -2WS[n],y = W^[n], and z = E x [n] - W^n]. 

The ISI model uses the values of W[n] and its deriva- 
tives at both the high-density (A — + 0) and the low- 
density (A — > 00) limits, to produce an accurate curve 



2 



for W(A), < A < 1, to insert in Eq. ([I]) to get an ap- 
proximation to E xc . It gives very accurate results for the 
correlation energy [7J and meets several conditions. But 
if we expand W lsl in the low density limit: 

W IS1 (X) = W 00 + ^=\-^ + ^A- 1 + • • • , (9) 
v Y i 

we can see that its A -1 term does not generally vanish, 
although it works very well numerically for E c |16j . This 
wrong coefficient was already shown to produce a spuri- 
ous term (A 2 In A) in the expansion of E c [n\] as A — ► oo 

There were several attempts to overcome this problem 
[correctly omitting the A~ x term but including all the 
other (integer and half-integer powers) terms] in the lit- 
erature [J2] [T7j by modifying the ISI model, but they 
are less simple: one requires Wq' [the next order in Eq. 

[17] and the other is not a direct model to W\ [T2] , 
Consider instead the following 4-parameter interpolation 
model: 

W acc (\)=a + by + dy\ y = —L=, (10) 

V 1 + cA 

where a,b,c, and d are scale-invariant functionals. We 
use the same inputs as those for the ISI model, i.e. 
Wo, Wo, Woo, an d W^,, to fit the parameters. Generally 
there are no analytical expressions in compact form for 
the parameters, and one has to solve for them numeri- 
cally. The 4th power in y is the lowest that can be added 
while satisfying the exact conditions, but producing an 
expansion with non-zero A - ™ terms (n € Z, n > 1). We 
recommend use of this W acc to replace the ISI model 
because it is numerically accurate and avoids the A -1 
term in the low-density limit. One can show that W" acc 
obeys the scaling property [Eq. Q], provided that 
W [n 7 ] = 7 W [n], WftnJ = Wftn], Wo> 7 ] = 7 Wx>H 
and W^o[n 7 ] = 7 3 / 2 W4,[n], as they should. If we in- 
tegrate W acc (A) over A from to 1, we find a simple 
expression for the exchange-correlation energy: 

Eg = a + ^ + 26(-l + VT+~c)/c. (11) 

We compare the performance of the new model and 
ISI on Hooke's atom, two electrons in a spherical har- 
monic well, with force constant k = 1/4. We show below 
that for this system, our W acc works as a highly-accurate 
interpolation, even more accurate than the ISI model. 

Magyar et al. [IB] calculated the W(A) curve for < 
A < 4 for Hooke's atom [k — 1/4) using Wq — E x — 
—0.515 and Wq = —0.101 as inputs. They confirmed 
that Woo = —0.743, consistent with the SCE ansatz [6]. 
They also found W^ = 0.235, but this was based on a 
fit that violated our condition, so we discount this result. 
Gori-Giorgi [19] calculated W'^ = 0.208 based on the 
SCE model [6] [12], which we consider exact. We apply 
these inputs (WhW^Wx., and W' x ) to our W acc and 
the ISI model {W acc generates two sets of solutions for 



a, b, c, and d, but we select the one with d closest to 

b, for it can be reduced to W slmp as below). We plot 
the differences between these models and the exact curve 
(taken from Ref. [TH]) in Fig. [I] One can see that our 
W acc works very well between A — and 1, which is the 
range of interest. Its predictions for W{,E C , and E c +T c 
are excellent, with T c being the correlation energy from 
the kinetic part, as listed in Table |TJ With these exact 
inputs, we found that, as A — ► oo, W lsl — > —0.743 + 

0.208A" 1 / 2 + 0.068A- 1 H , which shows that although 

the coefficient of A -1 is small, it does not vanish. 
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FIG. 1: Comparison of three different approximations to 
W{\) for Hooke's atom (Jfc = 1/4), plotted as AW = W model - 
IF cxact . The exact curve (up to A = 3) is taken from Ref. [18]. 



TABLE I: Comparison of several quantities for three different 
approximations to W(X) for Hooke's atom [k — 1/4). The 
exact values are taken from Ref. [18] except for W^ [19] . All 
energies are in mHartrees. 





exact 


ISI 


simp 


acc 


Wi 


-583 


-579 


-583 


-582 


w{ 


-44 


-41 


-45 


-44 


Ec 


-39 


-37 


-38 


-38 


Ec + Tc 


-10 


-10 


-9 


-9 



We can also apply our IF acc to the helium atom. Here 
Wq = E x = -1.025, W^ = -0.095 [2D], and W^ = 
-1.500 [6], Woo = °- 62 l [13 from the SCE model (SJUJ]. 
We plot the differences between these models and the 
exact curve (taken from Ref. [2T]) in Fig. [2] and compare 
several key quantities in Table |TT| 

One can see that our model here works fairly well, and 
W simp (see below) is even a little better than W acc . ISI 
does not satisfy the exact condition we derived in this 
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FIG. 2: Comparison of three different approximations to 
W(X) for helium atom, plotted as AW = W modcl - W cxact . 
The discrete values are shown, as well as fitting curves to aid 
the eyes. VV cxact values (up to A = f) are taken from Ref. 

ED- 



TABLE II: Comparison of several quantities for three different 
approximations to W(X) for helium atom. The exact values 
are taken from Ref. [21| . All energies are in mHartrecs. 





exact 


ISI 


simp 


acc 


Wi 


-1104 


-1100 


-1103 


-1103 


w[ 


-64 


-60 


-64 


-63 


E c 


-42 


-40 


-42 


-41 


E c +T c 


-6 


-6 


-5 


-5 



of W simp (A), and find: 
W - 



b = 



W„ 



4Wg 



(13) 



Thus a and b set the endpoints, while c is a measure of 
the curvature. Substituting Eq. (13 1 into Eq. (12 1, we 
get the explicit form of W(A) in terms of Wo, Woo and 
Wq. One can show that it has the correct expansion in 
both limits, and it obeys the scaling property [Eq. Q]. 
Setting d — b in Eq. (Ill and subtracting exchange, it 
yields: 



E s ™* = 2b[f(c)- 



f(c) = [VT+-c- 1 -±^]/c ) 
1 + c 



(14) 



with b and c defined in Eq. (13 1. £c" np correctly re- 
covers GL2 in the weakly-correlated limit (Woo — > — oo, 
keeping Wo and Wq fixed, such as in the Z — » oo limit 
of two-elecron ions) and E^ p correctly reduces to Woo 
for strong static correlation (Wq — » — oo, keeping Wo and 
Woo fixed, such as for stretched H 2 ). We can calculate 
the kinetic correlation energy T c : 



T c = b{2f(c)-z-z% 



(15) 



with /(c) defined in Eq. (14) and z = 1/yl + c, showing 
that the curvature j3 — T^J\E C — T c | 22J is a function of 
c alone. We strongly urge Exo be applied whenever its 
inputs are accurately known. 

We can further test our W slmp in systems with more 
than two electrons, but only those for which all inputs are 
known, with results listed in Table |III| One can see that 
W slmp predicts E c fairly accurately, but is less accurate 
than W ISI . This is perhaps due to lack of W' x in W simp . 



work [Eq. Q]: as A -> oo, W ISI -> -1.500+0. 621A~ 1 / 2 + 
0.376A -1 + • • • , so the A -1 coefficient is not even small. 

Now, we propose a simpler version of W acc , which can- 
not be used in typical cases, as the exact value of is 
not known in general. A simpler model is constructed by 
setting d = b, to yield: 

W sh ^(\) = a + b(y + y 4 ), y ^—L=^ (12) 

V 1 + cA 

with a, b and c being scale-invariant functionals. We have 
found (see results for Hooke's atom and helium atom) 
that although there is one parameter less, the above 
form produces usefully accurate results, especially be- 
tween A = and 1. In a word, W acc acts as an ac- 
curate interpolation to the whole adiabatic connection 
curve, while W slmp is more convenient and practical to 
use, without losing accuracy. It yields — 0.191 for 
Hooke's atom and 0.594 for helium. 

We use Wq, Woo and Wq to construct the explicit form 



TABLE III: Comparison of W simp and W ISI on systems with 
more than two electrons. Ex, Wo and Woo are taken from 
Ref. [13], and W^ is taken from Ref. \12\ . All energies are 
in Hartrees. 





Ex 




W SCE 

OO 


r^/SCE 

OO 




£Stmp 


rpcxact 

E c 


Be 


-2.67 


-0.250 


-4.02 


2.59 


-0.104 


-0.110 


-0.096 


Ne 


-12.1 


-0.938 


-20.0 


22.0 


-0.410 


-0.432 


-0.394 



In fact, in their first paper on the ISI model, Seidl et 
al. proposed a similar model [6], which yields results nu- 
merically very close to those of ISI, but without the y A 
term. But their model contains no \~ n (n > 1) contri- 
butions. Note that none of these models work for the 
uniform electron gas, because Wq = — oo [16], so both 
the model developed by Seidl et al. [B] and W slmp reduce 
to W(A) = Woo. 

After the bulk of this work was completed, we received 
a preprint of Ref. [12 , containing a detailed theory of 
the leading corrections to W(A) as A — > oo, consistent 
with the much simpler arguments given here. Also, we 
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use their value for helium (see text) to replace the We thank professor John Perdew for kind discussions, 
old one predicted by point-charge-plus-continuum (PC) This work is supported by National Science Foundation 
model [TJ. under grant number CHE-0809859. 
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